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Abstract 

We derive a new lower bound for the bandwidth of a graph that is based on a 
new lower bound for the minimum cut problem. Our new semidefinite programming 
relaxation of the minimum cut problem is obtained by strengthening the well-known 
semidefinite programming relaxation for the quadratic assignment problem by fixing 
two vertices in the graph; one on each side of the cut. This fixing results in several 
smaller sub-problems that need to be solved to obtain the new bound. In order to 
efficiently solve these subproblems we exploit symmetry in the data; that is, both 
symmetry in the min-cut problem and symmetry in the graphs. Our approach results 
in the best known lower bounds for the bandwidth of all graphs under consideration, 
i.e., Hamming graphs, 3-dimensional generalized Hamming graphs, Johnson graphs, 
and Kneser graphs, with up to 216 vertices. 

Keywords: bandwidth, minimum cut, semidefinite programming, Hamming graphs, John- 
son graphs, Kneser graphs 

1 Introduction 

For (undirected) graphs, the bandwidth problem (BP) is the problem of labeling the 
vertices of a given graph with distinct integers such that the maximum difference between 
the labels of adjacent vertices is minimal. Determining the bandwidth is NP-hard (see 
|35j ) and it remains NP-hard even if it is restricted to trees with maximum degree three 
(see [17] ) or to caterpillars with hair length three (see |34j ) . 

The bandwidth problem originated in the 1950s from sparse matrix computations, and 
received much attention since Harary's [20] description of the problem and Harper's paper 
[22| on the bandwidth of the hypercube. The bandwidth problem arises in many different 
engineering applications that try to achieve efficient storage and processing. It also plays 
a role in designing parallel computation networks, VLSI layout, constraint satisfaction 
problems, etc., see, e.g., [H|9l|32], and the references therein. Berger-Wolf and Reingold 
[2] showed that the problem of designing a code to minimize distortion in multi-channel 
transmission can be formulated as the bandwidth problem for the (generalized) Hamming 
graphs. 

The bandwidth problem has been solved for a few families of graphs having special 
properties. Among these are the path, the complete graph, the complete bipartite graph 
[6], the hypercube graph [23] , the grid graph [7], the complete k-level t-ary tree [39], the 
triangular graph [30], and the triangulated triangle [29]. Still, for many other interesting 
families of graphs, in particular the (generalized) Hamming graphs, the bandwidth is 
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unknown. Harper [24] and Berger-Wolf and Reingold [2] obtained general bounds for the 
Hamming graphs, but these bounds turn out to be very weak for specific examples, as our 
numerical results will show. 

We will consider the bandwidth of graphs with symmetry, such as the Hamming graphs, 
and use semidefinite programming (SDP) to obtain lower bounds, as it is known that such 
an approach is particularly fruitful in the case that the data have symmetry. Helmberg et 
al. |25j derived a lower bound for the bandwidth of a graph by exploiting spectral properties 
of the graph. The same lower bound was derived by Haemers |21j by exploiting interlacing 
of Laplacian eigenvalues. Povh and Rendl [37] showed that this eigenvalue bound can 
also be obtained by solving a semidefinite programming relaxation for the minimum cut 
problem. They further tightened the derived SDP relaxation and consequently obtained 
a stronger lower bound for the bandwidth. Blum et al. [3] proposed a SDP relaxation 
for the bandwidth problem that was further exploited by Dunagan and Vempalain [15] to 
derive an 0(log 3 n-y/log log n) approximation algorithm (where n is the number of vertices). 
De Klerk et al. [14] proposed two lower bounds for the graph bandwidth based on SDP 
relaxations of the quadratic assignment problem (QAP), and exploited symmetry in some 
of the considered graphs to solve these relaxations. Their numerical results show that 
both their bounds dominate the bound of Blum et al. [3], and that in most of the cases 
they are stronger than the bound by Povh and Rendl [37] . It is important to remark that 
all of the above mentioned SDP bounds are computationally very demanding already for 
relatively small graphs, that is, for graphs on about 30 vertices. 

In this paper we derive a new lower bound for the bandwidth problem of a graph that 
is based on a new SDP relaxation for the minimum cut problem. We prove that this new 
relaxation dominates the known SDP relaxations for the minimum cut. Consequently, we 
derive the best known lower bounds for the bandwidth of all graphs under consideration, 
i.e., Hamming graphs, 3-dimensional generalized Hamming graphs, Johnson graphs, and 
Kneser graphs. The computed lower bounds for Hamming graphs turn out to be stronger 
than the corresponding theoretical bounds by Berger-Wolf and Reingold [2], and Harper 
[24]. The new relaxation of the min-cut is a strengthened SDP relaxation of the QAP 
by Zhao et al. [43]. In particular, it is obtained from the latter relaxation by fixing two 
vertices in the graph; one on each side of the cut. Our approach results in several smaller 
subproblems that need to be solved in order to obtain the new bound. The number of 
subproblems depends on the number of 'types' of edges and non-edges of the graph, and 
this number is typically small for the graphs that we consider. In order to solve the SDP 
subproblems, we exploit the symmetry in the mentioned graphs and reduce the size of these 
problems significantly. We are therefore able to compute lower bounds for the bandwidth 
of graphs with as much as 216 vertices, in reasonable time. In order to evaluate the lower 
bounds, we also present upper bounds for the bandwidth of the above mentioned graphs 
by using a heuristic that improves the well-known reverse Cuthill-McKee algorithm [5]. 
Consequently, we are able to determine an optimal labeling (and hence the bandwidth) 
for several graphs under consideration. 

The further set-up of the paper is as follows. In Section [2] we introduce notation, 
provide some definitions, and give some background on symmetry in graphs. In Sections 
13.11 and 13. 2| we relate the bandwidth problem with the quadratic assignment problem and 
the minimum cut problem, respectively. In Section 13.21 we state Povh and Rendl's [37] 
expression for the lower bound for the bandwidth that exploits the lower bound for the 
minimum cut problem, and slightly improve this expression. The eigenvalue relaxation and 
the Povh- Rendl relaxation for the minimum cut are presented in Section [4J In Section [5J] 
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we exploit the fact that the minimum cut problem is a special case of the QAP, and suggest 
solving the well-known SDP relaxation for the QAP by Zhao et al. [13] to obtain a lower 
bound for the minimum cut. We also show that this SDP relaxation for the minimum 
cut dominates the relaxation by Povh and Rendl [37] . In Section 15.21 we derive a new 
SDP relaxation for the minimum cut problem that dominates all previously known SDP 
relaxations. In order to compute the new lower bound for the bandwidth of graph with 
symmetry, we exploit this symmetry as described in Section 15.31 Our numerical results 
are presented in Section [U where we also describe our heuristic for the upper bounds. 

2 Preliminaries 

In this section, we give some definitions, fix some notation and provide basic information 
on symmetry in graphs. 

2.1 Notation and definitions 

The space of k x k symmetric matrices is denoted by Sk and the space of k x k symmetric 
positive semidefinite matrices by S£ . We will sometimes also use the notation X y 
instead of J £ St, if the order of the matrix is clear from the context. For two matrices 
X, Y G R nxn , X > Y means > yij, for all i, j. The group of n x n permutation matrices 
is denoted by II n , whereas the group of permutations of {1,2,..., n} is denoted by Sym n . 

For index sets a, (3 C {1, . . . , n}, we denote the submatrix that contains the rows of A 
indexed by a and the columns indexed by (3 as A(a, f3). If a = /3, the principal submatrix 
A(a, a) of A is abbreviated as A(a). To denote column i of the matrix X we write X-^. 

We use I n to denote the identity matrix of order n, and to denote the i-th. standard 
basis vector. Similarly, J n and u n denote the n x n all-ones matrix and all-ones ra-vector, 
respectively. We will omit subscripts if the order is clear from the context. We set 
E%j — ^i^j ■ 

The 'diag' operator maps an n x n matrix to the n-vector given by its diagonal, while 
the 'vec' operator stacks the columns of a matrix. The adjoint operator of 'diag' we denote 
by 'Diag'. The trace operator is denoted by 'tr'. We will frequently use the property that 
tr AB = tiBA. 

For a graph G = (V,E) with \V\ = n vertices, a labeling of the vertices of G is a 
bijection (j) : V — > {1, . . . , n}. The bandwidth of the labeling (f> of G is defined as 

croo(G,0) := max \<j>(i) - <j>{j)\. 
0,j'}es 

The bandwidth a^G) of a graph G is the minimum of the bandwidth of a labeling of G 
over all labelings, i.e., 

o"oo(C) := min {<7oo(G, <fi) \ (j) : V — > {1, . . . , n} bijective} . 

The Kronecker product A(&B of matrices A S W xq and B £ W xs is defined as the pr x qs 
matrix composed of pq blocks of size r x s, with block ij given by aijB {i = 1, . . . ,p; j = 
l,...,q). The following properties of the Kronecker product will be used in the paper, 
see, e.g., [H] (we assume that the dimensions of the matrices appearing in these identities 
are such that all expressions are well-defined): 

(A®B) T = A T ®B T , (A® B)(C ® D) = AC ® BD. (1) 



3 



2.2 Symmetry in graphs 

An automorphism of a graph G = (V, E) is a bijection ir : 1/ — > V that preserves edges, that 
is, such that {tt(x), 7r(y)} G 22 if and only if {x, y} € E. The set of all automorphisms of G 
forms a group under composition; this is called the automorphism group of G. The orbits of 
the action of the automorphism group acting on V partition the vertex set V; two vertices 
are in the same orbit if and only if there is an automorphism mapping one to the other. 
The graph G is vertex-transitive if its automorphism group acts transitively on vertices, 
that is, if for every two vertices, there is an automorphism that maps one to the other (and 
so there is just one orbit of vertices). Similarly, G is edge-transitive if its automorphism 
group acts transitively on edges. In this paper, we identify the automorphism group 
of the graph with the automorphism group of its adjacency matrix. Therefore, if G 
has adjacency matrix A we will also refer to the automorphism group of the graph as 
aut(A) := {P G II n : P T AP = A}. 

As a generalization of the above algebraic symmetry, combinatorial symmetry is cap- 
tured in the concept of a coherent configuration as introduced by Higman [27] (see also 
[28J); indeed as a generalization of the orbitals (the orbits of the action on pairs) of a 
permutation group. It is defined as follows. 

Definition 1 (Coherent configuration). A set of zero-one n x n matrices {A\, . . . , A r } is 
called a coherent configuration of rank r if it satisfies the following properties: 

(i) X^iex = I for some index set X C {1, . . . , r} and Yll=i A% = J, 

(ii) Aj G {Ai, . . . , A r } for i = 1, . . . ,r, 

(iii) AiAj G span{j4i, . . . , A r } for all i, j. 

We call A := span{y4i, . . . ,A r } the associated coherent algebra, and note that this is 
a matrix *-algebra. If the coherent configuration is commutative, that is, A^Aj = AjAi 
for all i,j = 1, . . . ,r, then we call it a (commutative) association scheme. In this case, I 
contains only one index, and it is common to call this index (so Aq = 2), and d := r — 1 
the number of classes of the association scheme. 

One should think of the (nondiagonal) matrices Ai of a coherent configuration as 
the adjacency matrices of (possibly directed) graphs on n vertices. The diagonal matrices 
represent the different 'kinds' of vertices (so there are \I\ kinds of vertices; these generalize 
the orbits of vertices). In order to identify the combinatorial symmetry in a graph, one 
has to find a coherent configuration (preferably of smallest rank) such that the adjacency 
matrix of the graph is in A, see, e.g., [12]. In that case the nondiagonal matrices Ai 
represent the different 'kinds' of edges and nonedges. 

Every matrix *-algebra has a canonical block-diagonal structure. This is a consequence 
of the theorem by Wedderburn |42] that states that there is a *-isomorphism 

V : A — > ©f =1 C n > xn \ 

Note that in the case of an association schema, all matrices can be diagonalized simulta- 
neously, and the corresponding *-algebra has a canonical diagonal structure ©f =0 C. 

We next provide several examples of coherent configurations and association schemes 
that are used in the remainder of the paper. 
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Example 2 (The cut graph). Let m = (mi, 7712,1713) be such that mi + 7712 + = n. 
The adjacency matrix of the graph G miim2jrri3 , which we will call the cut graph, is given 
by 

(Omixmi J mi X m-j mi xm3 \ 
JmiV.m\ m2 xm2 m2 xm3 I • (2) 
0rrt3Xmi m3 xm2 m3 xm3 / 

The cut graph is edge-transitive, and belongs to the coherent algebra spanned by the 
coherent configuration of rank 12 that consists of the matrices 

J 0\ 

B 1 = | | , B 2 = | | , £ 3 = [ I = J3«f, 



/° J \ T 

5 4 = = fij, B 6 = 

\0 0/ 

B B = J J = Bj , B n = 
\0 0/ 

where the sizes of the blocks are the same as in ([2]). The coherent algebra is isomorphic 
to C © C © C © C 3x3 . In Appendix \A\ the interested reader can find how the associated 
*-isomorphism ip acts on the matrices Bi, for i = 1, . . . , 12. Note that the cut graph 
G mi mi m3 (that is, mi = m^) is in the coherent configuration of rank 7 consisting of the 
matrices B\ + B e , B 2 + B 7 , B 3 + B 5 , B 4 + B 8 , B 9 + B w , Bn, and Byi- 

Example 3 (The Hamming graph). The Hamming graph H (d, q) is the Cartesian product 
of d copies of the complete graph K q . The Hamming graph H(d, 2) is also known as the 
hypercube (graph) Qd- With vertices represented by d-tuples of letters from an alphabet 
of size q, the adjacency matrices of the corresponding association scheme are defined by 
the number of positions in which two (i-tuples differ. In particular, (Ai) XtV = 1 if x and y 
differ in i positions; and the Hamming graph has adjacency matrix A±. 

Example 4 (The generalized Hamming graph). The 3-dimensional generalized Hamming 
graph H qimm is the Cartesian product of K qi , K q2 , and K qs . With V = Q\ x Q 2 x 
Qz, where Qi is a set of size qi (i = 1,2,3), two triples are adjacent if they differ in 
precisely one position. Its adjacency matrix is the sum of three adjacency matrices in the 
corresponding 7-class association scheme on V, that is, the scheme with the 2 3 adjacency 
matrices describing being the same or different for each of the three coordinates. 

Example 5 (The Johnson and Kneser graph). Let fi be a fixed set of size v and let d 
be an integer such that 1 < d < v/2. The vertices of the Johnson scheme are the subsets 
of £1 with size d. The adjacency matrices of the association scheme are defined by the 
size of the intersection of these subsets, in particular (Ai) w w i = 1 if the subsets oj and u' 
intersect in d — i elements, for i = 0, . . . , d. The matrix A\ represents the Johnson graph 
J(v,d) (and Ai represents being at distance i in G). For d = 2, the graph G is strongly 
regular and also known triangular graph. 

The Kneser graph K(v, d) is the graph with adjacency matrix Ad, that is, two subsets 
are adjacent whenever they are disjoint. The Kneser graph K(5, 2) is the well-known 
Petersen graph. 
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3 Modeling the bandwidth problem 



3.1 The bandwidth and the quadratic assignment problem 

In terms of matrices, the bandwidth problem asks for a simultaneous permutation of the 
rows and columns of the adjacency matrix A of the graph G such that all nonzero entries 
are as close as possible to the main diagonal. Therefore, a 'natural' problem formulation 
is the following one. 

Let k be an integer such that 1 < k < n — 2, and let B = (bij) be the n x n matrix 
defined by 

{1 for \i — j\ > k 
' ' (3) 
otherwise. 

Then, if an optimal value of the quadratic assignment problem 

min tr(X T AX)B 
xen„ 



is zero, then the bandwidth of G is at most k. 

The idea of formulating the bandwidth problem as a QAP where B is defined as 
above was already suggested by Helmberg et al. and further exploited by De Klerk 
et al. [H]. In [13], two SDP-based bounds for the bandwidth problem were proposed; 
one related to the SDP relaxation for the QAP by Zhao et al. [43] . and the other one 
to the improved SDP relaxation for the QAP that is suitable for vertex-transitive graphs 
and that was derived by De Klerk and Sotirov [11]. The numerical experiments in [14] 
show that the latter lower bound for the bandwidth dominates other known bounds for 
all tested graphs. Although in [14] the symmetry in the graphs under consideration was 
exploited, the strongest suggested relaxation was hard to solve already for graphs with 32 
vertices. This is due to the fact that in general there is hardly any (algebraic) symmetry 
in B, that is, the automorphism group of B only has order two. 



3.2 The bandwidth and the minimum cut 

The bandwidth problem is related to the following graph partition problem. Let (Si, 52, £3) 
be a partition of V with \Si\ = rrii for i = 1, 2, 3. The minimum cut (MC) problem is: 

OPTmc := min J2 a ij 

ieS 1 ,j&S 2 

(MC) s .t. (S 1 ,S 2 ,S 3 ) partitions V 

\Si\ = nii, i = 1,2,3, 

where A = (a^) is the adjacency matrix of G. To avoid trivialities, we assume that mi > 1 
and m2 > 1- We remark that the min-cut problem is known to be NP-hard |18j . 

Clearly, if OPTmc > for some m = (mi, 7^2,^3) then a^G) > 1113 + 1, as noted in 
[3H 125] , Povh and Rendl |37] extended this by proving the following. 

Proposition 6. [37] Let G be an undirected and unweighted graph, and let m = (mi,m2,ms) 
be such that OPTmc > a > 0. Then 

cr(G) 00 > max{m 3 + l,m 3 + \V2a] - 1}. 
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Proof. Let <p be an optimal labeling of G, and let (Si, S2, S3) be a partition of V, such that 
4>(Si) = {1, . . . , mi} and 4>{S2) = {m\ + "13 + 1 5 • • • , n}. Let A be the maximal difference 
of labels over all edges connecting sets Si and S2 (so A < ^(G)) and let 5 = A — m 3 
(and note that 6 > 1). Now the result follows from the fact that the number of edges 
between Si and S2 is at most 5(5 + l)/2. Therefore, 5(5 + l)/2 > OPTmc> which implies 
that 5 > \y/2a] - 1. □ 

Proposition can be further improved by fully exploiting the above inequality on 5. We 
also first round a (see also [37]), which is possible because OPTmc is integer. Indeed, 
from 5(5 + 1) > 2 [a] and (Joo(G) > 771,3 + ^ one obtains the following corollary. 

Corollary 7. Xei G be an undirected and unweighted graph, and let m = (7711,7722,7723) be 
such that OPTmc > a > 0. Then 



o"oo(G) > m 3 + 



-W 2W+ i 



Thus, it follows that a strong bound on the bandwidth of a graph could be obtained from 
a good lower bound for OPTmc- Therefore, studying the min-cut problem is crucial for 
computing bounds for the BP. In the next sections we first present two known bounds for 
the min-cut problem, and then derive a new SDP relaxation for this problem. 



4 Eigenvalue bound and SDP bound for the minimum cut 

Helmberg et al. [25] derived a lower bound for the min-cut problem using the Laplacian 
eigenvalues of the graph. In particular, for m = (7774,777,2,777,3) this bound has closed form 
expression 

OPT ci g = --/i 2 A 2 - rMi^n, (4) 

where A2 and \ n denote the second smallest and the largest Laplacian eigenvalue of the 
graph, respectively, and \i\ and /J2 (with \i\ > 112) are given by 



- ^ — 77717772 ± \J Tn\Tn2(n — 777i)(?7 — 777,2)^ 



£H,2 

77 

Further, Helmberg et al. [25] concluded that if OPT e i g > for some m = (mi, 7772,7713) 
then a(G) > 7773 + 1. However, in this case also Corollary [7] may be applied to compute the 
bound for the BP. We remark that the bound on the bandwidth by Helmberg et al. [25j 
was also derived by Haemers |21| . 

Povh and Rendl [37] proved that OPT e ; g is the solution of a SDP relaxation of a certain 
copositive program. They further improved this SDP relaxation by adding nonnegativity 
constraints to the matrix variable. The resulting relaxation is as follows: 

min \ tr( D ® A)Y 

s.t. \ tr((Eij + Eji) ® I n )Y = rriiSij, l<i<j<3 

tr(J 3 Ea)Y = 1, 1 < i < n 

tr(V; T ® Wj)Y = mi, 1 < i < 3, 1 < j < n 

\ tr((Eij + Eji) <S> J n )Y = miinj, 1 < i < j < 3 
Y > 0, Y € S+ n , 



(MC C 



OPJ 
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where D = E 12 + E 21 G M 3x3 , V { = etuj £ M 3x3 , Wj = eju^ € R nxn , 1 < i < 3, 1 < j < n, 
and 5ij is the Kronecker delta. We use the abbreviation COP to emphasize that MCcop 
is obtained from a linear program over the cone of completely positive matrices. 

The numerical experiments in [37] show that MCcop can provide a significantly better 
bound for the bandwidth of a graph than the eigenvalue bound ([4]). We remark that in [37] 
the bound for the bandwidth of a graph based on MCcop is computed by using the result 
of Proposition [6l whereas for the eigenvalue bound the inequality (Toa{G) > + 1 was 
used. In Section [6j we compute the eigenvalue-based bound for the BP by using Corollary 

m 



5 A new bound for the minimum cut 

In this section we will present our new bound. In Section 15.11 we propose to solve the 
known SDP relaxation for the QAP by Zhao et al. [33] in order to obtain a lower bound 
for the MC. We will show in Theorem [9] that this relaxation dominates MCcop- I n Section 
15.21 we further strengthen it by fixing an arbitrary edge in the corresponding cut graph. 
This results in several subproblems that need to be solved to obtain the new lower bound 
for the min-cut problem, as stated in Theorem [TUJ In Corollary [15] we show that the 
strengthened bound dominates the bound of Section 15.11 In order to efficiently compute 
the new bound, we exploit symmetry in the data (that is, in the cut graph and the graphs 
for which we compute the min-cut/bandwidth) as described in Section [5.31 

5.1 The minimum cut and the quadratic assignment problem 

As noted already by Helmberg et al. [25], the min-cut problem is a special case of the 
following QAP: 

min -trX T AXB, (5) 

where A is the adjacency matrix of the graph G under consideration, and B (see ([2])) is 
the adjacency matrix of the cut graph G mijm2ima . Therefore, the following SDP relaxation 
of the QAP (see \A2>\ 136] ) is also a relaxation for the MC: 

min \tv(B ® A)Y 

s.t. tr(/ n <g> Ejj)Y = 1, ti(Ejj <g> I n )Y = 1, j = 1, . . . , n 
(MCqap) tr(/„ ® ( J n - I n )) + ( J n - I n ) (8) I n )Y = 

tr(jy) = n 2 

y>o, Yes\. 



ii- 



One may easily verify that MCqap is indeed a relaxation of the QAP by noting that 
Y := vec(X)vec(X) T is a feasible point of MCqap f° r A £ H n , and that the objective 
value of MCqap a t this point Y is precisely tr X T AXB. The constraints involving Ejj 
are generalizations of the assignment constraints Xu n — X u n — u n ; and the sparsity 
constraints, i.e., tr(/ n ® (J n — I n ) + (J n — /„) tg) I n )Y = generalize the orthogonality 
conditions AA T = A T A = /„. The constraint tr(jy) = n 2 ensures that there are n 
nonzero elements in the corresponding permutation matrix X. 

The matrix B, see ([2]), has automorphism group of order mi\m2\rn^\ when m\ ^ ni2 
and of order 2(mi\) 2 m^- when mi = va 2 - Since the automorphism group of B is large, 
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one can exploit the symmetry of B to reduce the size of MCqap significantly, see, e.g., 
|X2(, [TOl \TT\ 114] , Consequently, the SDP relaxation MCqap can be solved much more 
efficiently when B is defined as in ([2]) than when B is defined as in ([3]) (in general, the 
latter B has only one nontrivial automorphism). 

We will next show that MCqap dominates MCcop- I n order to do so, we use the 
following lemma from [36] that gives an explicit description of the feasible set of MCqap- 



Lemma 8. |36l Lemma 6] A matrix 
/ y(n) ... y(m) 



Y 



es+ y«) G R-x", iJ = 1> ... ?n> (6) 



i yO 1 ) . . . y(nn) 
is feasible for MCqap if an d only if Y satisfies 

(i) tr(J n (8 (J n - /„)) + {Jn - In) ® 4)^ = 0, 

(ii) try( w ) = 1 /or 1 < i < n, and £? =1 diag(F (ii) ) = u, 

(iii) u T F^') = diag(Y^')) T for 1 < i,j < n, and 

( iv ) Er=i yfe) = udiag(y( JJ ')) T /or 1 <j < n. 
Now we can prove the following theorem. 

Theorem 9. Let G be an undirected graph with n vertices and adjacency matrix A, and 
7711,7712,7713 > 0, mi + m >2 + 1^3 = n - Then the SDP relaxation MCqap dominates the 
SDP relaxation MCcop- 

Proof. Let Y 6 S + 2 be feasible for MCqap with block form (0). From Y we construct a 
feasible point Z £ S^ n for MCcop i n the following way. First, define blocks 

mi mi mi+m2 mi n 

= y Y^\ = J2 E Y&\ Z^ = Y^ ij \ 

i,j=l i=lj=mi + l t=l j=mi+m2+l 

mi+m,2 mi+iri2 n n 

Y y&\ z^ = y J2 y&\ z^ = y y^\ 

ij'=mi+l i=mi+l j'=mi+r?i2+l ij=mi+m2+l 




Z (U) Z (12) z (13) 
Z (21) Z (22) Z (23) 
Z (31) Z (32) Z {33) 

where ZW = (Z^) T for i < j. 

To prove that itr((£7jj + Eji) ® I n )Z = rriidij for 1 < i < j < 3, we distinguish the 
cases i = j and i ^ j. In the first case we have that 

mi 

tr(E« ® 7 n )Z = tr Z w = ^t r y (ii) = m h 

i=l 

where the last equality follows from Lemma [8] (ii). In the other case, we have that 
tr((£y + Sji) ® J n )Z = tr(Z^) + Z&>) = 0, 
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where the last equality follows from Lemma [8] (i) . 

The constraint tr(J3 ® Ea)Z = 1 for 1 < i < n follows from the constraint tr(I n ® 
EiijY = 1 for MCqap and the sparsity constraint. To show that tr(^ T <£> Wj)Z = for 
1 < i < 3, 1 < j < n, we will use Lemma[8] (iii) and tr(/ n (g) Ejj)Y = 1. In particular, let 
us assume without loss of generality that i = 1 and j = 2. Then 



tr(y^^ = W Wf UV 





From i^Y^u = 1 (see Lemma [8] (iii)) it follows that ti((Eij + Eji) C*D J n )Y = 2miirij for 
1 < i < j < 3. 

It remains to prove that Z y 0. Indeed, for every x £ M 3n , let x £ R n be defined by 

x T := <8) xj. n , u^ rl2 ® x Li u m 3 ® x J:n) ■ 

Then x T Zx = x T Yx > 0, since Y >z 0. Finally, it follows by direct verification that 
the objective values coincide for every pair of feasible solutions (Y, Z) that are related as 
described. □ 



Although our numerical experiments show that the relaxations MCqap an d MCcop 
provide the same bounds for all test instances, we could not prove that they are equivalent. 
We remark that we computed MCcop om y f° r graphs with at most 32 vertices, since the 
computations are very expensive for larger graphs. 



5.2 A new relaxation by fixing an edge 

In this section we strengthen the SDP relaxation MCqap by adding two constraints that 
correspond to fixing two entries 1 in the permutation matrix of the QAP ([5]). In other 
words, the additional constraints correspond to fixing an (arbitrary) edge in the cut graph, 
and an edge or a nonedge in the graph G. In order to determine which edge or nonedge 
in G should be fixed, we consider the action of the automorphism group of G on the set 
of ordered pairs of vertices. The orbits of this action are the so-called orbitals, and they 
represent the 'different' kinds of pairs of vertices; (ordered) edges, and (ordered) nonedges 
in G (see also Section T2.2p . Let us assume that there are t such orbitals Oh (h = 1,2, ... ,t) 
of edges and nonedges. We will show that in order to obtain a lower bound for the original 
problem, it suffices to compute t subproblems of smaller size. This works particularly well 
for highly symmetric graphs, because for such graphs t is relatively small. We formally 
state the above idea in the following theorem. 

Theorem 10. Let G be an undirected graph on n vertices, with adjacency matrix A, and 
t orbitals Oh (h = 1,2, ... ,t) of edges and nonedges. Let m = (mi, mi, rn^) be such that 
m i + m 2 + m 3 = n - Let {s\,S2) be an arbitrary edge in the cut graph G miim2jm3 (with 
adjacency matrix B as defined in ©j, and (r^i,?-/^) be an arbitrary pair of vertices in 
Oh (h = 1,2, . . . ,t). Let U n (h) be the set of matrices X £ Ii n such that X rhltSl = 1 and 
X rh2>S2 = 1 (h = 1,2, . . . ,t). Then 

mmtTX T AXB= min min trX T AXB. 
X£Tl n h=i,2,...,txen n (h) 
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Proof. Consider the QAP in its combinatorial formulation 

n 

From this formulation it is clear that if ir is an optimal permutation, then for every a € 
aut(^4) also air is optimal. Now let ir indeed be optimal, let h be such that (ir(si), vr(s2)) € 
Oh, and let <r be an automorphism of G that maps (ir(si),ir(s2)) to (r^i, r^)- Then 
7r* := (J7r is an optimal permutation that maps (si,s 2 ) to (r/j^r/^) and that has objective 

n 

^ a w * (i)^* = trX T AXB 

for a certain X £ II n (/i). This shows that 

min trX T ^4X.B > min min trXTAXB, 

xeu n h=i,2,...,txen n (h) 

and because the opposite inequality clearly holds, this shows the claimed result. □ 

Remark 11. In principle, the result of Theorem 1101 gives a lower bound for the min-cut 
of G for every edge (s\, s 2 ) in the cut graph G miim2im3 . However, this lower bound is the 
same for every edge because of the edge-transitivity of the cut graph. 

Remark 12. In [11], it is shown that in a QAP with automorphism group of A or B 
acting transitively (that is, at least one of the corresponding graphs is vertex-transitive), 
one obtains a global lower bound for the original problem by fixing one (arbitrary) entry 
1 in the permutation matrix X. 

In the following lemma, we show that when we fix two entries 1 in the permutation 
matrix X, we again obtain a QAP, but that is smaller than the original one. 

Lemma 13. Let X £ II n and ri,T2, s±,S2 £ {1, 2, . . . , n} be such that s\ 7^ S2, X ri81 = 1, 
and X r2jS2 = 1. Let a = {1, . . . , n}\{ri,r2} and (3 = {1, . . . , n}\{si, S2}, and let A and B 
be symmetric. Then 

tr X T AXB = tr X{a, /3) T {A(a)X(a, (3)B((3) +C{a,P)) + d, 

where 

C(a, P) = 2A(a, n)B( Sl ,P) + 2A(a, r 2 )B(s 2 ,(3), (7) 

and d — CL riri b SlSl -\- O r2 r2^S2S2 ^^•rxr^sxsi.' 

Proof. We will show this by using the combinatorial formulation of the QAP, splitting 
its summation appropriately (while using that A and B are symmetric matrices), and 
then switching back to the trace formulation (we will omit details). Indeed, let tt be the 
permutation that corresponds to X; in particular we have that vr(si) = r\ and tt{s2) = r 2 . 
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Then 



trX T AXB = (hr(i)w(j)bij 

a r 2 Tr(j)bs2j + ^ 

= tr X(a, /3) T A(a)X(a, fi)B(p) 

+ 2 tr X(a, /3) T (A(a, n)B(si, 0) + A(a, r 2 )B(s 2 ,P)) + d 

= ti X {a, /3) T (A(a)X (a, /3)B(/3) +d(a,/3))+d. □ 

Since A(a), B(/3) € S n ~2, where a, j3 are defined as in Lemma [T3"l and X(a,(3) £ n n _2, 
the reduced problem 

min trX T (A(a)X£(/3) + C(a,/3)) (8) 
xen n _2 

is also a quadratic assignment problem. Therefore computing the new lower bound for 
the min-cut problem using Theorem 1101 reduces to solving several SDP subproblems of the 
form 

H* h = min ±tr((£(/3) ® A{a h )) + Diag(c))Y + \d h 

s.t. tr(J (g) Ejj)Y = 1, tr(^- tg> I)Y = 1, j = 1, . . . , n - 2 
(MC£ X ) tr(/ <g> ( J - /)) + ( J - /) ® I)Y = 

tr(jy) = (n - 2) 2 
Y>0, YG5+_ 2)2 , 

where /i = l,...,t, I,J,Ejj G <S n -2, c = vec(C(a, /?)), /3 = {1, . . . , n}\{si, s 2 }, a' 1 = 
{1, . . . , n}\{rhi, ^2}) an d the constant d' 1 is defined as a! in Lemma [T3j Finally, the new 
lower bound for the min-cut problem is 

MCfi x = min at. 

h=l,...,t 

Thus, computing the new lower bound for the bandwidth problem MCfi x involves solving t 
subproblems. This seems like a computationally demanding and very restrictive approach, 
knowing that in general it is hard to solve MCg x already when n is about 15, see [38]. 
However, our aim here is to compute lower bounds on the bandwidth problem for graphs 
that are known to be highly symmetric. Consequently, for such graphs t is small and we 
can exploit the symmetry of graphs to reduce the size of MC^ X as described in Section IST31 
It remains here to show that the new bound MCfi x dominates all other mentioned 
lower bounds for the min-cut problem. To show this we need the following proposition. 
Note that we already know that the eigenvalue bound OPT c i g is dominated by MCcop 
(see [37] ) , which in turn is dominated by MCqap (see Theorem [9]) . 

Proposition 14. Let (si, s 2 ) be an arbitrary edge in the cut graph G mi>m2>m3 and (r^i, r^) 
be an arbitrary pair of vertices in Oh {h = 1, 2, . . . , t). Also, let a h = {1, . . . , n}\{r/ l i, r^} 
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and ft = {1, . . . , n}\{si, S2}- Then the semidefinite program 



min |tr(B ® ^4)y 



(MC^ AP ) 



s.t. tr(J n ® E^Oy = 1, tr(^j J n )Y = 1, 3 = 1,... 
tr(/ n (g) (J n - /„)) + (J n - I n ) ® I n )Y = 
tr(JF) = n 2 

tr(^ a( ® S rMirM )y = 1, i = 1,2 



y > o, y g 5 



is equivalent to MCg x (h = 1,2, ... , t) in £/te sense that there is a bijection between the 
feasible sets that preserves the objective function. 



Now, from Proposition O it follows indeed that the new bound dominates all others. 

Corollary 15. Let (si,S2) be an arbitrary edge in the cut graph G mim2m3 and (r/ l i,r/ l 2) 
be an arbitrary pair of vertices in Oh (h = 1,2, ... ,t). Then the SDP relaxation MCf} x 
dominates MCqap. 

In this section, we proved that the proposed new bound for the min-cut dominates the 
SDP bound MCcop, which in turn dominates OPT c j g . Numerical experiments by Povh 
and Rendl |37j and De Klerk et al. |14j show that it is already hard to solve MCcop f° r 
graphs whose size is 32. In the following section we show how one can exploit symmetry 
of the considered graphs to efficiently compute the new relaxation MCfi x . 

5.3 Reduction by using symmetry 

Computational experiments show that in general, the SDP relaxation of the QAP by 
Zhao et al. |43| . i.e., MCqap, cannot be solved in a straightforward way by interior point 
methods for instances where n is larger than 15, see, e.g., [38] . However, it is possible to 
solve larger problem instances when the data matrices have large automorphism groups, 
as described in \10\ 1 12 j . Further, in [11] \T3[ [T4"] it is shown how to reduce the size of the 
SDP relaxation that is obtained from the relaxation from |43j after fixing one entry 1 in 
the permutation matrix. Here, we go one step further and show how to exploit symmetry 
in the data to solve efficiently the SDP relaxation that is obtained from MCqap after 
fixing two entries 1 in the permutation matrix, i.e., MCg x . 

In |1CH [T2] it is shown that if one (or both) of the data matrices belong to a matrix 
*-algebra, then one can exploit the structure of the algebra to reduce the size of the SDP 
relaxation by Zhao et al. [43] , In particular, since the matrix B belongs to the coherent 



Proof. The proof is similar to the proof of Theorem 27.3 in [40] , 



□ 
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algebra described in Example [2j MCqap reduces to the following: 

min ±trA(X 3 + X 5 ) 
s.t. X 1 +X 6 + X 11 = I n 
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Xi — j n 



i=l 



tr(JXi)= Pi , Xi>0, * = 1, . . . , 12 



l ^ 0, X 6 i-rl 7 >- 0, X u l -^X 12 y 

1 mi — 1 z — ' ° ni2- 1 ' — ' 11 m3 — 1 lz — 



/■±-(Xi + x 2 ) -A — x 3 



(9) 



l —x A \ 

urns * 



^3 = ^51 ^4 = -X9 > -^8 = -Xioi 
-^2, -XV > -^11, ^12 £ <->ra, 

where pi = mi, j>2 = mi (mi - 1), Ps = P5 = mim 2 , P4 = P9 = Wim3, y>6 = rn 2 , 
P7 = m 2 (m 2 - 1), p 8 = Pio = rri2m3, pn = m 3 , and pi 2 = m 3 (m 3 - 1). To derive ©, we 
applied the associated *-isomorphism given in Appendix [Al 

It is possible to perform further reduction of ([9]) for graphs with symmetry, see, e.g., 
In order to do so, we need a coherent algebra containing A. For the graphs that 
we consider, such coherent algebras are known, see Examples [SHU Therefore, further 
reduction of ([9]) follows directly from results from, e.g., [TOl [TTj. 

The symmetry reduction of the SDP relaxation MCg x (h = 1, . . . , t) is however not a 
straightforward application of results from [10 ^ I1 H [13 | I14j . Therefore we provide a detailed 
analysis below. In order to perform the desired symmetry reduction of MCg x , we will find 
a large enough subgroup of the automorphism group of the data matrix from the objective 
function, i.e., of B((3) ® A(a h ) + Diag(vec(C')), where C is given in ((7]). Note that if the 
objective function would not have a linear term, then aut (-£?(/?)) ® aut(^4(a ft- )) could serve 
as such a subgroup. 

In what follows, we first determine appropriate subgroups of the automorphism groups 
of A(a), B(f3), and C(a,/3). To describe these we introduce the following definition. 

Definition 16. For 7*i,r 2 G {1, . . . , n} with r\ 7^ r 2 , the subgroup of aut(A) that fixes 
row and column r\ and row and column r 2 of A is: 

stab((n,r 2 ), A) := {P G aut(A) : P ri>ri = P r ^ r2 = 1}. (10) 

This is the pointwise stabilizer subgroup of aut(A) with respect to (ri,r 2 ), see, e.g., [4]. 

It is important to distinguish the pointwise stabilizer from the setwise stabilizer of a 
set S C V; in the latter it is not required that every point of S is fixed, but that the set 
S is fixed (that is, n(S) = S). For a = {1, . . . , n}\{ri, r 2 } we define 

H(A(a)) := {P(a) : P G stab((n, r 2 ), A)}. 

This group is a subgroup of the automorphism group of A{a). 

Similarly, we define H(B{p)) := {P{fi) : P G stab((«i, s 2 ), B)} for p = {1, . . . , n}\ 
{s±, s 2 }. In the following, we will consider an edge {s\, s 2 } in the cut graph (without loss of 



h 
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generality s± G Si and s 2 G 5*2). Because of the simple structure of the cut graph, 7i(B(/3)) 
can be easily described. Indeed, because si and s 2 are fixed by P 6 stab((si, s 2 ), B), the 
sets 5i \ {si} and 52 \ {s 2 } are fixed (as sets) by P(/3) G H(B((3)). This implies that 
H(B((3)) is the direct product of the symmetric groups on Si \ {si}, S 2 \ { s 2}, and £3. 
In fact, this is the full automorphism group of B(/3) in case mi 7^ m 2 (in the case that 
mi = 1^2 it is an index 2 subgroup of aut(f?(/3)) since the 'swapping' of Si \ {si} and 
S2 \ { s 2} is not allowed). In any case, the action of %(P(/3)) has 12 orbitals (as descibed 
in Example [2]) . 

Now we can describe the group that we will exploit to reduce the size of the subproblem 
MC£ X (h = l,...,t). 

Proposition 17. Let n,r 2 , si,s 2 G {1, . . . , n}, a = {1, . . . ,n}\{n,r 2 }, /3 = {1, . . . ,n}\ 
{si,s 2 }, andC(a,f3) = 2A(a,r 1 )B(s 1 ,(3) + 2A(a,r 2 )B(s 2 ,f3). Then 

H(B(p))®H{A(a)) 
is a subgroup of the automorphism group of B(/3) <8> A(a) + Diag(vec(C(a, (5))). 

Proof. Let Pg G %(B{(3)) and Pa G %(A{a)). It is clear that Pg ® P4 is an automor- 
phism of P(/3) <8> A(a), so we may restrict to showing that is also an automorphism of 
Diag(vec(C(a, /?))). In order to show this, we will use that for i = 1,2 we have that 

PjTHag(A(a,ri))P A = Diag(A(a, n)) and p£Diag(B(/3, Sl ))P B = Diag(P(/3, Si )). 

Indeed, the first equation is equivalent to the (valid) property that a^uy. = aj ri for all 
j ^ ri,r 2 and all automorphisms n of A that fix both ri and r 2 , and the second equation 
is similar. Because vec(C(a, /?)) = 2B(j3,s{) ®A(a,ri) + 2B(/3, s 2 ) ®A(a,r 2 ), the result 
now follows from 

(P B ® P A ) T Diag [B(0, si) ® A(a, n) + P(/3, s 2 ) A(a, r 2 )} (P B ® P A ) 
= pTDiag(P(/3, Si ))Pb PjDiag(^(a, n))P A 

+ pTDiag(P(/3, s 2 ))Pb ® PjBmg(A(a, r 2 ))P A 
= Diag[P(/3, Sl ) ® A(a, n) + P(/3, s 2 ) A(a, r 2 )], 
where we have used the properties in ([1]) of the Kronecker product. □ 

Now, for an associated matrix *-algebra of A{a) one can take the centralizer ring (or 
commutant) oiH(A(a)), i.e., 

A A[a) = {X G R nxn : XP = PX, VP G n(A(a)}. 

Similarly, we take the centralizer ring A-b(8) of %{B (/?)). Now we restrict the variable 
Y in MCg x to lie in A.B(p) ® Aa{<x) > anci obtain a basis of this algebra from the orbitals 
of H(_A{a)) and 7i(B(f3)). Finally, the symmetry reduction of MCg x is similar to the 
symmetry reduction of MCqap, see ([9]) and also [TT| [T3], 

6 Numerical results 

In this section we present the numerical results for the bandwidth problem for several 
graphs. All relaxations were solved with SeDuMi [51] using the Yalmip interface [33] on 
an Intel Xeon X5680, 3.33 GHz dual-core processor with 32 GB memory. To compute 
orbitals, we used GAP [16]. 
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6.1 Upper bounds 

In order to get more information about the quality of the lower bounds, we needed good 
upper bounds for the bandwidth of the presented examples. We obtained these by using 
MATLAB's implementation of the reverse Cuthill-McKee algorithm [5|. However, the 
output seemed far from optimal, so we improved the results in the following way. 

Consider a labeling of the graph with bandwidth Coo. A vertex u is called critical if 
it has a neighbor w such that \4>{u) — 4>(w)\ = a^. Now we consider the critical vertex u 
that has the largest label, and its critical neighbor w. Let z be the vertex with the largest 
label that is not adjacent to any of the vertices with labels 1,2,..., 4>{w). If 4>(z) < 4>(u), 
then we decrease by one the labels of the vertices with labels <p(z) + 1, 4>{z) + 2, . . . , <p(u), 
and give z new label 4>(u). It is not hard to show that this operation does not increase the 
bandwidth of the labeling. We keep repeating this (to each new labeling) until we find no 
z such that cf)(z) < (f)(u). 

Because the Cuthill-McKee algorithm requires an original labeling (its input is the 
adjacency matrix of the graph), and the output depends on this labeling, we first permuted 
the labeling randomly, and repeated the entire algorithm 1000 times. The obtained output 
seems to be of good quality: in many cases for which we know the optimal value, this value 
is attained; for example for the Hamming graph H(3, 6), where we have an upper bound 
101, which was shown to be optimal by Balogh et al. [T]. We remark that the best obtained 
upper bound for H(3, 6) computed by MATLAB's Cuthill-McKee algorithm after 1000 
random starts was 130, while the above described improvement reports the optimal value 
(101) 22 times (out of 1000). Also for the Johnson graphs J(v, 2) with v G {4, . . . , 15} 
the improved Cuthill-McKee heuristic provides sharp bounds. Although the required 
computational times are acceptable (computing the upper bound for H(3, 6) takes 50 
seconds, for example), we do not provide details for all the graphs, because we did not 
do any attempt to optimize the implementation. We also do not claim that the above is 
a good heuristic that can compete with the state of the art. The purpose was to have 
a quick and dirty way to obtain reasonably good upper bounds in order to compare our 
lower bounds. 

6.2 Lower bounds 

In this section we present several lower bounds on the bandwidth of a graph. Each such 
lower bound, i.e., bw c ; g , bwcoP) bwQAP , and bwfi x is obtained from a lower bound of 
the corresponding relaxation (OPT e i g , MCcop> MCqap, and MCfi x , respectively) of some 
min-cut problem. Indeed, for each graph and each relaxation we consider several min-cut 
problems, each corresponding to a different m. In particular, we first computed OPT c i g 
and the corresponding bound for the bandwidth (using Corollary [7]) for all choices of 
m = (mi, 7712,7713) with mi < 777-2; this can be done in a few seconds for each graph. The 
reported bw e i g is the best bound obtained in this way, and m e i g is the corresponding m. 
We then computed MCqap (or MCcop f° r the hypercube graph) for all m with mi < mi 
and 7773 > 7773 (where m = m e i g ). Similarly we computed the bound on the bandwidth 
that is obtained from MCg x . 

The hypercube graph. The first numerical results, which we present in TablefU concern 
the Hamming graph H(d,2), also known as the hypercube Qd, see Example [3j Because 
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the bandwidth of the hypercube graph was determined already by Harper [23] as 

this provides a good first test of the quality of our new bound. In the other families of 
graphs that we tested, we did not know the bandwidth beforehand, so there the numerical 
results are really new (as far as we know). Besides that, the tested hypercube graphs are 
relatively small, so that we could also compute MCcop (which we do not do in the other 
examples) . 

Table [1] reads as follows. The second column contains the number of vertices n = 2 d of 
the graph, and the last column contains the exact values for the bandwidth (Joo{Qd)- in 
the third column we give the lower bound on the bandwidth of the graph corresponding 
to OPT c ig, while in the fourth-sixth column the lower bounds correspond to the solutions 
of the optimization problems MCcop, MCqap, and MCg x , respectively (obtained as de- 
scribed above). We remark that when the lower bound bwQAp was tight (this happened 
for d = 2,3), we did not compute the bound bwg x , since then this is also tight. Note that 
for d = 2,3, 4, the latter bound is indeed tight. 



d 


n 


bw eig 


bw C op 


bwQAp 


bw fix 


Poo(Qd) 


2 


4 


2 


2 


2 


2 


2 


3 


8 


3 


4 


4 


4 


4 


4 


16 


5 


6 


6 


7 


7 


5 


32 


7 


10 


10 


11 


13 



Table 1: Bounds on the bandwidth of hypercubes Qd- 



In Table O we list the computational times required for solving the optimization prob- 
lems MCcop, MCqap, and MCfi x that provide the best bound for the bandwidth. The 
computational time to obtain MCfi x for Qd is equal to the sum of the computational times 
required for solving each of the subproblems MCg x , h = 1, . . . , d. Table [2] provides also the 
best choice of the vector m, i.e., the one that provides the best bandwidth bound for the 
given optimization problem. We remark that for d = 3, 4 there are also other such best 
options for m. Table [11] (see Appendix [B]) gives the number of orbitals in the stabilizer 
subgroups H(Qd(&)), for different a = {1, . . . , n}\{ri,r2}. 



d 


?7i e ig 


MCcop 


mcoP 


MCqap 


"IQAP 


MC fix 


w fix 


2 


[1,2,1] 


0.51 


[1,2,1] 


0.05 


[1,2,1] 






3 


[3,3,2] 


1.29 


[2,3,3] 


0.53 


[2,3,3] 






4 


[6,7,3] 


23.55 


[7,4,5] 


1.32 


[7,4,5] 


7.12 


[4,6,6] 


5 


[11,15,6] 


1019.83 


[10, 14,8] 


1.86 


[10,14,8] 


74.56 


[10, 12, 10] 



Table 2: Qd- time (s) to solve relaxations and corresponding m. 



The Hamming graph. Next, we give lower and upper bounds on the bandwidth 
of the Hamming graph H{d, q) for q > 2, see Example [3l Because the two-dimensional 
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Hamming graph H(2, q) (also known as the lattice graph) has bandwidth equal to ^ q+ 2 ' q — 1 
(see |26j). we computed bounds on the bandwidth for the next two interesting groups of 
Hamming graphs, i.e., H(3,q) and H(4,q). Tables [3] El and 1121 are set-up similarly as in 
the case of the hypercube graphs, except for the last column of Table [3] where instead of 
the exact value of the bandwidth we now give the upper bound obtained by the heuristic 
as described in Section 16.11 

Concerning this upper bound, we remark that it matches the exact bandwidth for 
H(3,6) (see pQ). Besides this case, the upper bounds for the Hamming graph from the 
literature are weaker than ours. For example, Harper's upper bound \24\ Claim 1] for 
H(3, 4) is 33, while the upper bound for H (3, 5) by Berger-Wolf and Reingold [2j Thm. 2] 
is 84. 

For the lower bounds, it is interesting to note that bw e i g equals bwQAP for all cases 
H(3, q) with q = 3, . . . , 6, while the new lower bound is better for all graphs in the table. 

Concerning Table H] we note that for all instances there are more options for m that 
provide the best bound. The number of orbitals in , H([H(4,3)](a)) for each subproblem 
as given in Table [12] is large, but we are able to compute MCfi x because the adjacency 
matrix A(a) has order only 79. 



d 


q 


n 


bw eig 


bwQAP 


bw fix 


u.b. 


3 


3 


27 


10 


10 


12 


13 


3 


4 


64 


22 


22 


25 


31 


3 


5 


125 


43 


43 


47 


60 


3 


6 


216 


74 


74 


78 


101 


4 


3 


81 


21 


23 


26 


35 



Table 3: Bounds on the bandwidth of the Hamming graphs H(3,q) and H(4,q). 



d 


Q 


TTleig 


MCqap 


"1QAP 


MC fix 




3 


3 


[9,10,8] 


0.29 


[9,10,8] 


44.03 


[4,12,11] 


3 


4 


[21,22,21] 


2.80 


[21,22,21] 


176.38 


[20, 20, 24] 


3 


5 


[42, 43, 40] 


15.50 


[42, 43, 40] 


536.65 


[38,41,46] 


3 


6 


[72, 74, 70] 


76.20 


[72, 74, 70] 


1756.88 


[63, 76, 77] 


4 


3 


[27,35,19] 


9.10 


[29,30,22] 


5877.33 


[22, 34, 25] 



Table 4: H(3,q) and H(4,q): times (s) to solve relaxations and corresponding m. 

The 3-dimensional generalized Hamming graph. In Tables [5] [6j and [13] we 

present the analogous numerical results for the 3-dimensional generalized Hamming graph 
H qi) g 2t q 3 , as defined in Example H] To the best of our knowledge there are no other lower 
and/or upper bounds for the bandwidth of H gimm in the literature. Our results show 
that the new bound can be significantly better that the eigenvalue bound. For instance, 
the eigenvalue lower bound on the bandwidth of #3,4,5 is 16, while bwfi x is 24. 

The Johnson graph. In TablesEKU and 1 141 we present the analogous numerical results 
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<?1 


Q2 


93 


n 


bw e ig 


bwQAp 


bw fix 


u.b. 


o 

I 


3 


3 


18 


5 


8 


9 


9 


2 


3 


4 


24 


6 


10 


11 


12 


2 


3 


5 


on 

3U 


1 


ii 


13 


15 


2 
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Table 5: Bounds on the bandwidth of H, 
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Table 6: H qi 

mm'- time (s) to solve relaxations and corresponding m. 



for the Johnson graph J(v,d), as defined in Example [5) In particular, we provide bounds 
for J(v, 3), with v € {6, . . . , 11}, and J(8, 4). The bandwidth of the Johnson graph J(v, 2) 
(also known as the triangular graph) has been determined by Hwang and Lagarias |30j . 
and equals L^ 2 /4J + \v/2] — 2. We remark that when the bound bwQAP was tight (this 
happened for v = 6,7), we did not compute the bound bwg x , since then the latter is also 
tight. Indeed, it thus follows that the bandwidth of J(6,3) equals 13 and the bandwidth 
of J(7, 3) equals 22. 

We also remark here that there are more options for m that provide the best bound, 
in particular for J(u,3), with v > 8. For example, the lower bound for the bandwidth 
of J(8,3) is equal to 31 for the vectors [12,14,30] and [13,13,30]. Table O see Ap- 
pendix [HI provides the number of orbitals from the stabilizer subgroups Tl{[J(v, d)](a)), 
a = {1, . . . , n}\{ri,r2} for d = 3 and d = 4, respectively. Since for d = 4 the number of 
orbitals increases significantly when v increases from eight to nine, we could not compute 
the new lower bound bwfi x ( J(9, 4)). However, we obtained that bw e i g ( J(9, 4)) = 49, and 
bw QA p(J(9,4)) = 52 in 23.12 seconds. 
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Table 7: Bounds on the bandwidth of J{v , 3) and J(v,4). 
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[32,32,56] 


11 


3 


[43, 57, 65] 
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1607.29 


[38,56,71] 
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[21,22,27] 


2.31 


[13,26,31] 
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Table 8: J(v,3) and J(t>,4): time (s) to solve relaxations and corresponding m. 

The Kneser graph. In Tables [9] and [TO] we finally present our bounds and compu- 
tational times for the Kneser graph K(v, 2) with v € {5, . . . , 8} and K(v, 3) with v £ 
{7, . . . , 10}, see Example [5j We remark that the orbitals of &ut(K(v,d)) and aut(J(u,d)) 
are the same, and the corresponding orbitals in the stabilizer subgroups for these two 
graphs are also the same. 

We note that Juvan and Mohar [31] presented general lower and upper bounds for the 
Kneser graph. Their lower bound however is an eigenvalue bound that is weaker than the 
eigenvalue bound bw e i g by Helmberg et al. [25] . Also their general upper bound is (much) 
weaker than our computational results. 

Since the bound bwQAP for the Petersen graph K(5, 2) is tight, we did not compute 
bwfi x for this graph. Note that it is a folklore result that the bandwidth of the Petersen 
graph equals 5. 
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A *-isomorphism from Example [2] 
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Orbitals in stabilizer subgroups 
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Table 11: Number of orbitals in l~L{Qd{ a ))- 
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Table 12: Number of orbitals in H([H(d,q)](a)). 
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Table 13: Number of orbitals in T~L{H qitq2) q 3 (a)) 
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Table 14: Number of orbitals in U([J(v, 3)] (a)) and W([J(u,4)]( 
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